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' We consider the feasibility of using the Galileo Navigation Satellite System to con- 

^\ , strain possible extensions or modifications to General Relativity, by assessing the impact 

, of the related additions to the Newtonian potential and comparing with the available 

' observables: the relative frequency shift and the time delay of light propagation. We 

address the impact of deviations from General Relativity based on the parameterized 
Post-Newtonian parameters due to the presence of a Cosmological Constant, of a con- 
stant acceleration like the putative Pioneer anomaly, a Yukawa potential term due to 
massive scalar fields and a power-law potential term, which can arise from Ungravity or 
f(R) theories. 
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1. Introduction 

The Galileo positioning system is an important step towards the improvement and 
development of new applications in navigation monitoring and related topics. Its 
operational use of precision clocks in orbit and comparison with those on ground 
stations enables one to view it broadly as a timing experiment in outer space. Hence, 
Galileo offers a great opportunity for fundamental research in physics: together with 
the deployed Global Positioning System (GPS) and GLONASS systems, satellite 
navigation is indeed the first technological application where relativistic effects are 
taken into account as a regular engineering constraint on the overall design (see 
Refs. ^'^'"^ and references therein). 

As such, there are several effects arising from special and General Relativity 
(GR) that must be taken into account, i.e. time dilation, gravitational blueshift and 
the Sagnac effect. These may yield a clock deviation of as much as ^ 40 fis/day: 
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many orders of magnitude above the precision of the onboard clocks considered in 
the mentioned navigations systems, which is of the order of 4 ns/day, i.e. a time 
stability of about 5 x lO^^'^. 

Furthermore, the gravitational frequency shift is of the order of Vn/c^ — 10"-'^° 
(where Vn = GMe/Re is the Newtonian potential, G is Newton's constant. Me — 
6.0 X 10^^ kg is the Earth's mass. Re — 6.4 x 10^ m is its radius and c is the speed 
of light). In the Galileo Navigation Satellite System, this correction is accounted by 
the receiver, while the GPS system accounts for this mismatch through an offset in 
the onboard clock frequency. 

The Galileo navigation system offers a positioning improvement of at least one 
order of magnitude, from an everyday error margin of ~ 10 m with the GPS system 
to ~ 1 to; a spatial accuracy of the order of 1 toto is possible using the planned thirty 
satellite constellation and carrier phase measurements, similar to the increased pre- 
cision of the GPS system using averaging and modern geodetic GPS receivers. Given 
this, a legitimate question arises: what are the possible implications for fundamental 
physics that one may extract from this increased precision? 

In this study, we aim to establish bounds on the detectability of extensions and 
modifications to GR, by assessing the impact of the related additions to the Newto- 
nian potential on the observables made available by the Galileo system, namely the 
relative frequency shift e/ = ///q — 1 and the light propagation time delay At. This 
may be done for a wide variety of already available phenomenological models (see 
Rcfs. for updated surveys). Hence, this study should be regarded as an extension 
of the excellent treatments already available that deal with physical effects arising 
from GR i'^-^-^, but in the context of extensions of the latter. 

We assume as typical values Cf = 10"-'^^ (a precision possible using next gen- 
eration Rubidium clocks) and A.t = 10~^^ s (corresponding to a 1 toto precision); 
throughout the text, results stemming from the comparison with At = 10~^ s will 
also be given, so the reader may grasp the observable range obtained for spatial 
accuracies ranging between 1 m and 1 toto. 

The use of the global positioning system to ascertain the effects induced by 
modifications of gravity here assessed would involve the direct comparison between 
the received signal and its emitted time and frequency. However, one rcicalls that 
the actual frequency emitted by the satellites is not known, but may be inferred 
from comparison with clocks at the ground segment, assuming that they behave 
similarly (which appears reasonable, given the great stability of the onboard and 
ground clocks). 

In this approach, one does not consider the increased precision due to averag- 
ing procedures or interpolation of the signals transmitted by the full constellation 
of satellites. A more evolved approach could resort to the rcidundancy introduced 
by the various elements of the constellation to eliminate systematic errors, correct 
for atmospheric and propagation effects, validate the comparison with the ground 
clocks, etc.. Regarding the latter point, one remarks that by comparing several 
received signals with the master clocks on the ground, one could lift the degener- 
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acy between "true" frequency shifts and those arising from a slow, unpredictable 
drift away from the base rate of the emitter clocks (constrained by their ~ 10~^^ 
stability). 

A more evolved use of the multitude of available satellites would rely on the 
construction of the so-called GPS (or emission) coordinates, whereas one resorts to 
four time coordinates, instead of one time plus three spatial ones, to define a suitable 
frame of reference '''•^■^•lo. These are obtained via four material bodies exchanging 
light signals, and are thus naturally implemented by global positioning systems. 
Prom a theoretical point of view, they may be regarded as a realistic implementation 
of Einstein's point coincidence timing procedure. 

The realization of these GPS coordinates requires solely that one is able to 
measure the direction and arrival time of the emitted signal, as well as its redshift 
— precisely the quantities we will use as observables in this study. In this fashion, 
the GPS coordinates allow one to directly read the metric components in the related 
frame and, via a suitable transformation, to another frame {i.e. one expressed in 
terms of time as measured by an observer and the usual spherical coordinates). 

After this, one could then compare with the metric (in a rotating frame) arising 
from the approximately spherical mass distribution of the Earth, and interpret 
the differences as due to fundamental modifications of gravity. The accuracy and 
stability of the on-board clocks is reflected in the uncertainty of this transformation: 
again, combinations of the available satellites into sets of four could be used to 
obtain a group of GPS coordinates and observables that might help to eliminate 
systematics and other sources of error. 

This paper is organized as follows: firstly, we briefly review the main relativistic 
effects that arc accounted for in the Galileo navigation system. We then proceed 
and consider the possibility of measuring several corrections to the law of gravity 
using Galileo: 

• Deviation from GR based on the parameterized Post- Newtonian parameters 

• Presence of a Cosmological Constant 

• Constant acceleration like the putative Pioneer anomaly 

• Yukawa addition mediated by massive scalar fields 

• Power-law addition, which can arise from Ungravity or f{R) theories 

Finally, conclusions are drawn and an outlook is presented. 
2. Main relativistic effects 

As stated, the purpose of this study is to assess the feasibility of the Galileo system as 
a probe of extensions of GR. In order to better contextualize this, the effects arising 
from this theory that are relevant for the implementation of a global navigation 
satellite system are briefly reviewed, according to discussions that can be found 
in Refs. i'2'7.3 -pjjg excellent discussion of Ref. 1 on GR and the GPS system is 
followed, with the suitable adaptations and computations for the Galileo system. 
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for comparison. The reader who is familiarized with this discussion may skip these 
introductory paragraphs and proceed directly to the new results found from Section 
3 onward. 



2.1. Frame of reference 

One begins by assuming that all time-dependent effects are of cosmological origin, 
and thus evolve over a timescale of order Hq^, where Hq is Hubble's constant; 
hence, one may discard these as too small within the timeframe of interest and 
assume a static, spherically symmetric scenario. Given this, one considers the stan- 
dard solution of the Einstein field equations for a non-spherically symmetric mass 
distribution — which, in isotropic form, is given by the line element 



ds'= 1 



2V 



(c dtf 



1 



2V 



{dr^ + dfl^) 



(1) 



= (l + ^)(cdif-(l-^)(dr^+dfl^) , 

where dQ^ = r'^{d9'^ +sm^ OdcjP) is the solid angle element, and V is the gravitational 
potential. In the unmodified GR scenario, this is simply the Newtonian potential 



V = VN{r,e,ct>) = - 



GMe 
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where 9 and (f) arc the latitude and longitude and Pnm arc Legendrc polynomials of 
degree n and order m; J„ are the zonal harmonic coefhcients, independent of the 
longitude; for n 7^ to, the quantities C„m and Snm are the tesseral harmonic coef- 
ficients, while for n = TO they are dubbed sectoral harmonic coefficients. The three 
types of harmonic reflect the mass distribution of the Earth and its deviation from 
sphericity For very low eccentricity orbits, such as those found in global position- 
ing systems, the tesseral and sectoral harmonics give rise to periodic perturbations 
only 

One now takes into account the rotation of the Earth with respect to this fixed- 
axis reference frame, with angular velocity cj = 7.29 x 10~^ rad/s; the so-called 
Langevin metric may be obtained by performing a coordinate shift t' = t, r' = r, 
6' = 9 and ^' = <p — cot', yielding the line element 



ds' = 



i9 



2wr^ sin^ 9d(l)dt 
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(dr 
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where, for simplicity, primes were dropped. 

Clearly, a non-diagonal element appears, plus an addition to the gravitational 
potential interpreted as a centrifugal contribution due to the rotation of the refer- 
ence frame; this leads to a definition of an effective potential ^ = V — (ursinO)'^ /2. 
The parameterization of the Earth's geoid is obtained by taking the multipole ex- 
pansion of V up to the desired order, calculating the equipotential lines $ = $o 
(the value of $ at the Equator) and solving for r{9,(f)). The value adopted by the 
International Astronomical Union is $o/c^ = -6.969290134 x 10~^°. 

In the above, the coordinate time t is equal to the proper time of an observer at 
infinity. Given the issue of ground to orbit clock synchronization, the metric should 
be rewritten with a time coordinate coincident with the proper time of clocks at 
rest on the Earth's surface. 

Since the already discussed geoid provides one with an equipotential surface 
= all clocks at rest with respect to it beat at the same rate. Hence, rescaling 
the time coordinate as t {1 + $o/c^)i yields 



1 + 



2($ - $c 



(c dtf 



(4) 



2ujr^sin^ed(l)dt -{^- ^^"^ {dr^ + d^^) ■ 

In the above, a similar transformation of the spatial coordinates (realized in the 
usual scheme of GPS and Galileo) is not performed, for simplicity. Going back to 
an inertial, non-rotating frame, one finally writes the metric as 



ds' = 



1 + 



2(1/ -$o) 



(c dtf 



2V 



(5) 



2.2. Constant and periodic clock deviation 

In this paragraph one discusses the difference in clock rates due to relativistic effects, 
i.e. the difference between coordinate time t as measured in an Earth centered, non- 
rotating frame, and the proper time r of the moving clocks onboard the satellites. 

Keeping only terms of order c~^, the proper time increment on the moving clock 
is approximately given by 



dr = ds/c=(l+^^^^^^dt , (6) 

so that, assuming an elliptic orbit with semi-major axis a and, for simplicity, the 
Newtonian potential generated by a perfectly spherical body V = Vn — —GMe/v, 
one obtains ^ 



dr = ds/c- 



3GMe _ £o _^ 2GMe 
2a(? (? (? 



dt 



(7) 
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Hence, the constant correction terms are given by 

SGMe $0 in 

+ ^ = -4.7454 X 10-10 , (8) 

for Galileo, and —4. 4647 x 10^^", for the GPS. The difference stems from the slightly 
higher orbit for the former, hcaiueo ~ 23.2 x 10^ km against hcps ~ 20.2 x 10^ km. 
As the above shows, the orbiting clock beats faster by about 41 us/day (Galileo) 
and 39 fis/day (GPS). 

Other residual periodic corrections arc proportional to (1/r — 1/a); one may 
directly integrate the corresponding term in Eq. (7), yielding ^ 

^tperiodic = -^vGMEasinE = 4.4428 x 10~'^°ey/asmE sj^/m , (9) 

where E is the eccentric anomaly, defined by the transcendental equation E — 
esini? = ^GMeI a?{t — tp), with tp the time of perigee passage. Assuming an 
almost circular orbit with eccentricity e = 0.01, this contribution can reach 23 ns 
(GPS) or 24 ns/day (Galileo, due to the slightly higher orbit) — see Ref. 13 for a 
discussion on the long-term deviation from circular orbits. 

These second-order effects arise from Eq. (7), via the substitution of the New- 
tonian expression 

into Eq. (7); since the purpose of this study is to ascertain the effect of additional 
corrections to the Newtonian potential Vjv, it is important to determine the effect 
of these additions in this computation. 

To do so, one considers a spherically symmetric addition U{r) to the Newtonian 
potential, V = Vn + U. When a satellite is orbiting at r = a, the derived force is 
centripetal; together with the assumption energy conservation (since the additional 
contribution is assumed to be conservative), this yields 

^ = K(a) + C/» , (11) 

^ + VN{r) + U{r) = ^ + VN{a) + U{a) . 

Considering for simplicity the spherically symmetric Newtonian potential Vjv = 
—GMe/t, this leads to 



with Ua = U{a) and f/„ = U'{a), for brevity. 



(12) 
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Before replacing Eq. (12) into Eq. (6), one should notice that the additional term 
U will also shift the value of the effective potential at the Equator, $o ^ + Ue 
(with Ue = U{ai), whore oi is the equatorial radius). Considering this, one obtains 
the generalization of Eq. (7), 



dr = ds/c = 



1- 



3GMe 3>g 



1 fU'a 



2ac2 

2GMe fl 

a 



+ UE-Ua] + 



(13) 



[U{r) - Ua] 



dt 



Solving for dt and taking only c ^ terms, one obtains 



At 



I dr 

J path 



3GMe 
2ac2 



r.2 



2 



Ue - Ua 



(14) 



2GMe 



U{r)] 



Hence, perturbing the Newtonian potential yields an additional constant rate cor- 
rection [U^a/2 + Ue — Ua\/(?, as well as a periodic term 2{Ua — U)/(?. 



2.3. Doppler Effect 

The transmission of a light signal between two moving bodies affects the frequency 
fs of the GPS signal by the usual special relativistic (first-order) Doppler effect, as 
well as the already mentioned second-order effect in Eq. (6). Since the correction 
mechanism only takes the constant clock rate shift {2^o+3GME)/2ac'^ into account, 
the received frequency /e will be shifted with respect to the corrected frequency 
fs by the linear Doppler effect, the unaccounted additional constant term and the 
Newtonian and perturbative periodic terms. 



/e = fs 



^ Ka ^ Ug-UE ^ 2GMe fl 1 
2c2 



2(C/ - Ua) 



n • u 



n ■ V 



(15) 



where v is the velocity of the satellite, u = Re x lu the rotation velocity of the 
geoid and n a unit vector in the direction of propagation (as seen in a local inertial 
frame). 

As seen in Eq. (12), if an addition to the Newtonian potential is regarded as a 
perturbation, then the velocity of a satellite orbiting at a radius r is given by 



IGMe f2 



Ua-U- 



2 



GMp 



(16) 



Following Fig. (1), the deviation of the Doppler factor at the end of Eq. (15) is 
given by 
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Fig. 1. Signal path geometry: tlic gray area (not to scale) gives the ~ 1 km? area where the 
linear Doppler term is negligible, 6 <g 1. "Nose down", y = 0, and grazing transmissions, y ~ Re, 
are also depicted (dotted and full line, respectively). 



= , „ , - ^ 3 X 10-Q-f- . (17 

c + f^y c Rs + hc Re 

The maximum value is obtained for a grazing signal y ^ Re (since h > Re)- the 
perturbation U will of course modify its magnitude, due to the shift in the velocity. 

Given this, comparison of Eq. (16) with (15) indicates that the linear Doppler 
shift affects the contributions of the constant clock rate term —U'^a/2c^ and the 
periodic term [U — Ug)/ ■ At leading order, the relative frequency shift induced 
by U, dubbed e/, is 



2c2 c- 



,2 



2{u - Ua) 



2 

with the term arising from the linear Doppler shift being defined as 



(18) 



r/ ^ _ y URE + h)c^ 

^^y^=R^\l GMe ^''^ 

setting a = Re + h. 

Clearly, the competition between the second-order effect embodied in the metric 
Eq. (14) and the linear effect of the kinematic Doppler factor depends upon the 
angle between the velocity of a satellite and the direction of propagation, i.e. the 
ratio u/Re- one may ascertain the values for which the linear effect cancel the 
second-order one in Eq. (18): 
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and S{y) = 2 ^ ?y - 10"'' i?^. Hence, one concludes that only for an almost radial 
transmission ( "nose down" , with a spread of less than 1 km) is the linear term 
smaller than the second-order one. 

One thus assumes the scenario of a grazing transmission y = Re. where the linear 
term is enhanced and the possibility of detecting an addition to the gravitational 
potential is boosted. A more realistic choice would involve a shorter transmission 
path y < Re, thus avoiding the issue of transversing the thick atmospheric layers 
and deal with the ensuing signal distortion. Also, the relation bc;twcen the different 
orbital planes adopted by each global navigation satellite constellations and the 
rotation axis of the Earth would have to be properly accounted for, namely through 
a careful evaluation of the Sagnac effect for each case. 

With the above in mind, one writes 



ef = ec + ep , (21) 



Prom Eq. (20), one sees that the parameter S{Re) decreases with y/h; since 
hoaiiieo ~ 23.2 X 10^ km against hcps ~ 20.2 x 10^ km, one has 6{RE)Gaiiieo = 
1.7 X 10"* versus S{Re)gps = 1-8 x 10"*. Thus, the GPS global navigation system 
yields a marginally higher frequency shift (a 5% increase of the affected terms), 
although insufficient against the argued advantage of using the European system 
increased precision to detect it. 

Finally, notice that the constant frequency shift Cc in Eq. (21) differs significantly 
from the naive expectation that would arise if one straightforwardly considered the 
Schwarzschild metric and simply computed the gravitational frequency shift due to 
U (setting Vjv = 0), 



I goo Earth I - 2U E / , Ua - Ue 

The difference is not a minor one: the dominance of the linear Doppler term implies 
that the induced frequency shift can be several orders of magnitude higher than the 
above would anticipate. 

In the subsequent sections, one shall use the characteristic values for this con- 
stellation, with the notation 5 = S{RE)Gaiiieo, for brevity; also, one approximates 
the equatorial radius ai by the mean Earth radius Re, and the semi-major axis as 
a = Re + h, and computes the amplitude of the periodic clock rate term, given by 
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Cp = 6{Re) 



Ug-Ub 
C2 



(23) 



where Ub = U{b) and b = ay/1 — is the semi-minor axis. 
2.4. Shapiro time delay and the Sagnac effect 

The Shapiro time delay is a second order relativistic effect affecting the propagation 
of light ^, given by 



a result obtained after integrating over a radial path of proper length / and con- 
sidering only the effect of the Newtonian potential. In the following sections the 
additional contribution stemming from the specific perturbative potentials U will 
be addressed; for simplicity, one uses I = h, thus assuming a "nose-down" propaga- 
tion path. 

The Sagnac effect is yet another relativistic contribution, which reflects the rota- 
tion of the Earth and the consequent difference between the gravitational potential 
y of a non-rotational frame and the effective potential $ of a rotating one; from 
the standpoint of an observer in an inertial {i.e. non-rotating) frame, it reflects the 
added propagation delay due to the rotation of the receiver as light travels towards 
it. From Eq. (4), one gets the additional time delay 



where the integral gives the equatorial projection of the area swept out by the 
light signal in its travel between the satellite transmitter and the rotating ground 
receiver. Its magnitude it highly dependent on the assumed path: given the typical 
configuration of ground segment stations and the location of the satellites, it varies 
between 240 ns and 350 ns . 

As already shown above, a global positioning system is affected by a frequency 
shift of order 10^^" and a cumulative propagation time delay of the order of 0.1 [is 
(mainly due to the Sagnac effect). The following sections aim to compute other 
corrections that should be taken into consideration when calculating the clock syn- 
chronization of ground and onboard clocks, and compare the obtained results with 
the precision of 10~^^ and the time accuracy of Galileo, of order 10~^^ s. 

3. Post-Newtonian effects 

Before venturing into more speculative and hypothetical effects, arising from ex- 
tensions to GR, it is natural to first tackle the possibility of measuring Post- 
Newtonian effects with the Galileo positioning system; these are naturally much 




(24) 




(25) 
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smaller than the previously considered, of higher order in the Newtonian potential, 
Vn/c^ !v GMe / (Rec^) ~ 10"-^° (again considering a purely spherical body). 

Post-Newtonian effects are suitably addressed by resorting to the so-called pa- 
rameterized Post-Newtonian (PPN) formalism, which allows one to describe higher- 
order effects induced by metric extensions and alternatives to GR. For simplicity, 
one focuses only on the (3 and 7 PPN parameters, thus writing the related PPN 
metric 5,i5,i6,i7 



ds^ = - 



2V (V 



(26) 



The /3 parameter measures the amount of non-linearity affecting the superposition 
law for gravity, while 7 is related to the spatial curvature per unit mass. 

Although not evoked here, the full PPN metric includes a total of ten PPN 
parameters, which characterize the underlying fundamental theory and its possible 
consequences; these may include a violation of momentum conservation, existence of 
a privileged reference frame, amongst others deviations from GR. The PPN formal- 
ism is defined so that General Relativity is parameterized by = — 1, while all 
remaining parameters vanish; measurements of the Nordtvedt effect yield |/3 — 1| < 
(2 — 3) X 10"'', while Gassini radiometry indicates that 7 — 1 = (2.1 ± 2.3) x 10~^ 

18,19 

Also, one recalls note that the gravitational potential V appearing above is, in 
principle, due to the matter and energy distribution of the whole Universe, through 



V{x) 



^ I p{x',t) 
J \x-x'\ 



d^x' 



(27) 



with similar definitions for the metric potentials involved in the full PPN met- 
ric (with the aforementioned ten parameters), which include not only the density 
distribution, but also the velocity field of the source matter (see Ref. 5 for a full 
discussion). 

Notice that the PPN potentials do not depend on the PPN parameters, only 
on the matter density, pressure and velocity: the latter are coefficients affecting the 
expansion of the metric elements in terms of these quantities, as Eq. (26) shows. One 
may neglect the contributions of distant sources to these potentials, independently 
of possible deviation of the PPN parameters from their GR values (experimentally 
constrained to be very small); in other words, no PPN parameter directly enforces a 
Machian view of gravitation, in which distant sources would directly induce strong 
local effects on the gravitational field. 

With the above in mind, one writes the acceleration in the weak-field limit, 



^00 ~ ~ 29^^900 



GMp 



2(7-^) 



(28) 
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where the prime denotes differentiation with respect to the radial coordinate. 

The experimental bounds discussed above indicate that the difference (7 — /?) 
should be of order < 10^'*; since GMe/c?R ~ 10^^", one concludes that PPN 
corrections to the acceleration arc much too small to be detected. Likewise, sim- 
ilar expressions may be derived for the time delay and frequency shift, showing 
that the post-Newtonian relative corrections are indeed proportional to Vn/(? = 
GMe/{Re<?) ^ 10^^". Given the already considered gravitational time delay and 
frequency shift are of the order Aj ~ 10~^ s and e/ ~ 10"^", respectively, compar- 
ison with the available precision of Galileo and GPS systems makes it clear that 
Post-Newtonian effects signaling deviations from GR are much below the observa- 
tion threshold. 

3.1. Strong Equivalence Principle 

There is another impact of deviations from GR within the PPN formalism that must 
be considered. Indeed, while distant bodies do not impose a strong acceleration on 
the individual satellites, the effect of the small gradient between their effect on the 
latter, and on the Earth itself — that is, a tidal acceleration — is boosted if the 
Strong Equivalence Principle (SEP) is not valid This difference would be perceived 
on the ground as an acceleration of the satellite with respect to the geocentric frame, 
and could in principle be relevant. Thus, the purpose of this section is to ensure 
that such an effect is under control in the context of the approach of this work. 

The SEP implies that the gravitational field properties refiect the gravitational 
energy of the bodies themselves, and provides an assumption about the nonlinear 
properties of gravity. GR assumes that this principle is exact, but alternative met- 
ric theories of gravity, such as those involving scalar fields, and other extensions. 



In the framework of the PPN formalism (and assuming fully conservative, 
Lorentz-invariant theories, for simplicity ^'^''), the validity of the SEP translates 
into the vanishing combination ry = 4/3 — 7— 3 = — trivially satisfied for GR, 
where /3 = 7 = 1. If 7^ 0, the acceleration of a body with inertial and gravitational 
mass m and nig, respectively, is given in the Newtonian limit by 



typically violate the SEP 20,21,22,23 




(29) 



The ratio between gravitational and inertial masses is given by 



rUg 



E 



(30) 



1 



7] 



.2 



TO 



mc- 



with the (negative of) the gravitational self-energy E given by 




(31) 
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Let US now consider two bodies with inertial mass mi, m2 in the gravitational 
field of an external object with mass M at a distance L; the bodies arc separated 
by a small distance I when compared to L. For simplicity, one considers only the 
case where the three bodies are aligned, whereas the tidal acceleration, given by the 
difference between individual gravitational pulls, attains its higher value: 



02 - fli = GM 
GM " 



1 



mg2 



m2 



L2 



(32) 



i2 

GM 
1? 



mi 

21 7] 



21 



mi {L-iy 

m2 _ 
E2 El 
^m2 mi 
using I <^ L. 

The term utn = 2GMl/L^ is the Newtonian tidal acceleration. Considering the 
Earth and one satellite, mi = Me and m2 = m « 700 kg, separated by a distance 
I — Re + /i « 3 X 10^ m, one may compute the ratio between this term and the 
gravitational pull of the Earth on a satellite, ue = GMe/I"^- 



ClTN 
Ue 



Me 



(33) 



One can obtain this ratio for the relevant bodies in the Solar System: the Sun (M = 

3.3 X 10^ Me and L = 1 AU) yields aTN/o-E ^ 10~^; Jupiter at its closest point 
from Earth (M = 318Me and L = 4.2 AU) produces arN/aE ~ 10"^°; finally, the 
Moon (M = 1.23 x IQ-^Me and L = 3.84 x 10^ m) leads to axN/aE ~ 10"^ Thus, 
one concludes that the tidal forces induced by external bodies may be disregarded 
if the SEP holds and r/ = 0. 

This said, there is the possibility that SEP violating effects might induce addi- 
tional, relevant tidal effects, as expressed by the second term of Eq. (32), 



am = V 



GM 



Ee 



Es 



(34) 



where Ee = 4.6 x 10~ ^"M^jc^ is the gravitational self-energy of the Earth ^, and 
Es the corresponding quantity for a navigation satellite. Indeed, one notes that the 
Newtonian tidal acceleration drops for increasingly distant external bodies (as L^), 
while the term arising from SEP violation falls only quadratically, thus eventually 
surpassing the former. 

The extreme smallness of gravitational self-energy for laboratory-sized objects'* 
implies one has to consider astrophysical scales in order to test the SEP. Currently, 



''It is of the same order of magnitude as a sphere with radius R ^ I m and uniform density, 
yielding Es/mc^ = (3/5) (Gm/iic^ ) ~ IQ-^s. 
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the Earth — Moon — Sun system provides the best scenario for testing this principle, 
manifested in the so-called Nordtvedt effect Lunar laser ranging experiments 
yield the constraint -i] = (4.4 ± 4.5) x 10^^ 

One can easily check that no Solar System body can lead to a SEP-violation 
induced term larger than the Newtonian tidal acceleration qtn- Indeed, the ratio 
between the two is 



^«^^^<2xl0-^3^: . (35) 

ttTN MeC^21 I ^ ' 

Taking the highest allowed value yields the upper bound ij < 8.9 x 10^'*. 

The SEP violation induced term is important only if the above ratio is much 
larger than unity, i.e. when the distance between the external body and the Earth- 
satcllitc system is L ^ 5 x 10^^/ ~ 5 kpd Thus, one concludes that only objects at 
a galactic distance might yield relevant effects. 

The above result only states that for galactic ranges, the SEP violating tidal 
term is larger than the Newtonian one, aTri > citn, but not necessarily measurable. 
To assess this, one considers the ratio between the former and the gravitational pull 
of the Earth as, 

^«^^f}V#.<4xlO-"^fl)' ' (36) 
as Me \L J Mec^ Me \L J ^ ' 

again considering the highest allowed value for r]. Considering the Milky Way mass 
to be concentrated at its core, M = Mmw = 2.3 x lO'^^Af^ and L = 7.62 kpc, one 
obtains aTrj/dE < 1-6 x 10~^^; Virgo, the closest cluster to Earth, at a distance 
L = 16.5 Mpc and with a mass M = 4.0 x 10^°Mb, yields ar^/aE < 5.7 x 10~^^. 

One now recaps the results of this section: although external bodies at a suffi- 
ciently large distance (at the galactic range) may induce SEP violating tidal effects 
larger than their Newtonian counterpart, these are clearly much below the grav- 
itational pull of the Earth: one finds that the external bodies are not sufficiently 
massive to compensate for their extreme remoteness. As a result, one can safely 
state that SEP violating effects can be disregarded in the context of the present 
study; conversely, one does not have to assume that this fundamental tenet of GR 
is obeyed, i.e. the discussion presented before concerning the PPN parameters is 
not constrained by the condition 7/ = 4^ — 7 — 3 = 0. 



4. Detection of the cosmological constant 

According to the latest observations, the Universe is currently ongoing a period of 
accelerated expansion; although several proposals exist to account for this acceler- 
ation, the simplest explanation resorts to a cosmological constant A 10~^^ s~^, 
which acts as a fluid with negative pressure (see Ref. 26 and references therein). The 
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local effect of this component may be evaluated by matching the outer Friedmann- 
Robertson- Walker metric with a static, symmetric solution given by Birkhoff 's the- 
orem; this yields the Schwarzschild-de Sitter metric^, with a line element ^8,29^ 

-1 

in anisotropic form; the presence of a cosmological constant produces an additional 
potential term Ua = -Ar^/6 (with A = SQaH^ and Qa = 0.7). 

In order to compare with the metric Eq. (5), a coordinate change to an isotropic, 
co-rotating frame of reference should be performed. However, given the weak field 
regime, this would only amount to a small correction to the effect of this additional 
potential, which is in itself a perturbation to the Newtonian one. 

4.1. Constant frequency shift 

Following Eq. (21), the constant frequency shift induced by this extra potential 
contribution is easily obtained, 

ec~-(5A(i?^ + /j)2^_10-34 , (38) 
clearly much below the e/ = 10~^^ precision available. 

4.2. Periodic frequency shift 

The maximum periodic frequency shift is given by 

Cp = -S-^{Re + hfe^ = CcC^ < ec < e/ , (39) 
and so is impossible to detect. 

4.3. Propagation time delay 

Similarly, the additional propagational time delay is given by 

Atdelay = ^ / UA{r) dr , (40) 
c JRe 

so that the presence of a cosmological constant results in a further contribution of 



ds^ = - 1 + 



2{Vn + Ua) 



(c dtf + 



1 + 



2{Vn + Ua) 



^As a side note, one remarks that the identification of the time coordinate t of the current form 
as the proper time of an observer at rest at infinity breaks down, due to the collapse of the 
Schwarzschild "bubble" at a finite distance, where it matches the exterior FRW metric 
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1 Ar^ A 

= ^ y -^dr= [{3Re{Re + h) + h''] ~ lO'^^ s , (41) 

also many orders of magnitude below the optimum time resolution of 10~^^ s. 
Therefore, one concludes that the cosmological constant is completely undetectable 
by the Galileo positioning system (as indicated by an analytical study in Ref. 30). 



5. Detection of anomalous, constant acceleration 

Although not usually considered, the presence of an anomalous, constant, accelera- 
tion affecting the free-fall of bodies could model effects arising from some fundamen- 
tal theory of gravity, perhaps hinting at the existence of a fundamental threshold 
between known dynamics and yet undetected exotic physics. 

One widely discussed example is the so-called Modified Newtonian Dynamics 
(MOND) model 3i>32,33^ which features a departure from the classical Poisson equa- 
tion at a characteristic acceleration scale of about 10^^" m/s^, and aims to explain 
the puzzle of the galaxy rotation curves without evoking any dark matter compo- 
nent. 

From the experimental viewpoint, a constant acceleration ap = (8.74 ± 1.33) x 

10"-'^'' m/s^ is reported to affect the Pioneer 10/11 probes 34:35,36. j-j-g Qj-jgjj^^ either 
due to an incomplete engineering analysis (sec e.g. Ref. 37, 38, 39, 40, 41, 42) or 
stemming from yet undiscovered fundamental physics 43,44,45^ j^g^g |-,ggji dubbed the 
Pioneer anomaly. 



5.1. Constant frequency shift 

An anomalous, constant acceleration, ac, would imply an addition to the gravi- 
tational potential of the form Uc = acf] following the procedure depicted in the 
previous section, one calculates the related constant frequency shift as 



ecc^5^{RE + h) , (42) 

and, comparing with the frequency stability e/ = 10~^^, one obtains an upper bound 
for a detectable constant acceleration ac, 



6 {Re + h) 

This range allows for a narrow fit of the reported Pioneer anomalous acceler- 
ation ap, which yields ec = 2.5 x 10~^^; notice that previous computations 
without the linear Doppler shift yielded a completely negligible ec, which shows the 
enhancing effect of the former. 
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Since a more evolved treatment is not pursued here (including realistic locations 
of ground segment stations, atmospheric correction, etc.) and this unobserved cor- 
rection is at the edge of detectability, one cannot regard this finding as a major 
point against of the existence of the Pioneer anomaly — which, given the conclu- 
sions of Refs. 37, 38, 39, 40, 41, 42, now presents an unaccounted acceleration of 
about only one third its originally reported value ap. 

5.2. Periodic frequency shift 

Similarly to the result of paragraph 4.2, the maximum periodic frequency shift is 
given by 

2 

ep~(5^^(i?£; + /i) = ece^<ec , (44) 

and is therefore of no use in further constraining the allowed range of the anomalous 
acceleration ac- 

5.3. Propagation time delay 

The propagational time delay due to this extra potential addition is given by 

^tc=- ^dr= ^h{2RE + h) . (45) 

By the same token, comparison with a time accuracy of 10~^^ s indicates that only 
accelerations a > 0.1 m/s^ are measurable using time delay. 

6. Detection of Yukawa potential 

A Yukawa potential is one of the more discussed modifications to the law of gravity, 
as it may arise from scalar-vector-tensor field "fifth force" models, where its char- 
acteristic range A is related to the mass m of the scalar or vector field, A a 

4 

For the case of exchange of a scalar particle, the full potential is given by 

V{r) = -^^(l + ae-^'>') , (46) 

where a < 1 is the strength of the perturbation and Gqo the gravitational coupling 
for r oo; the latter redefines Newton's constant G through the relation G = 
Goo(1 + q;); this full potential may be separated into a Newtonian-like potential and 
an extra potential 



(47) 
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By conjugating several constraints arising from different setups covering a wide 
range of distances (from near-millimeter tests to planetary experiments), stringent 
bounds have been obtained for the allowed region of parameter space a, A, as 
may be seen in Fig. 2. 

10' 
10^ 
10-3 
10* 
10 = 

m 10* 

10-' 
10* 

lO-to 

10^2 10^ W 10* 10« l& 10» 10** w* 

Fig. 2. Exclusion plot (shaded region) for a Yukawa-type additional force with strength a and 
range A: dark gray corresponds to the available bounds *® and light gray gives the superimposed 
excluded range obtained for a relative frequency precision = 10~^^; approximations in the 
short- and long-range regimes are depicted as gray lines. 

One now computes the constant relative frequency shift. 

Since the magnitude of A is not known, terms in S may not dominate this expression. 
To ascertain this, one must consider the short- and long-range regimes separately. 
The max;imum periodic term is given by 




J I I \ I I I I I I I I I I I I L 



(Re + h)c^ 
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(49) 

6.1. Short-range "fifth force" 
6.1.1. Constant frequency shift 

One first considers that the additional Yukawa interaction is short-ranged, A ^ 
h, Re- In the regime A < Re, the Yukawa potential term Eq. (47) is further 
affected by a form factor that reflects the extended mass distribution of the source 
(the Earth) and the cixchange of massive particles of mass m = A~^ only within a 
spherical shell of limited width. 

This form factor serves to further suppress the already exponentially small efl:ect: 
the constant frequency shift Eq. (48) is approximated by 

ReC^ 

and the Doppler linear term is negligible for 



. uGMe 
'jR^Thj^' 



-(RE+h)/X 



(50) 



2Ae''/^ > SRe X< OMRe , (51) 

Thus, one does not pursue the computation of the effects of a short-range Yukawa 
potential. Fig. (2) depicts the numerical computation of the constraint arising from 
the constant frequency shift with no form factor included: one concludes that no 
bound on the parameter space (A, a) can be obtained in this regime, even with this 
unrealistic assumption. 



6.2. Long-range "fifth force" 
6.2.1. Constant frequency shift 

One now examines the opposite assertion concerning the characteristic Icngthscalc 
A, and instead assumes a long range fifth force, A ^ h,RE- The exponential terms 
may be expanded to first order in r/A and the additional frequency shift Eq. (48) 
is dominated by the linear Doppler shift term'^, 



.aGMERs + h 
^c--^^^^ A^ • (^2) 

Notice that the zeroth-ordcr term arising from the expansion of the potential Eq. 
(47) is absorbed into the constant G. 



°As noted before, a previous computation 46,47 -^yithout this term yielded a value for ec smaller 
by a factor i5 ~ 10'', again displaying the relevance of this effect. 



July 19, 2011 0:17 WSPC/INSTRUCTION FILE GalileoIJMPD4 



20 O. Bertolami and J. Paramos 



For this effect to pass unnoticed at a e/ = 10 ° precision, one must have 
log" « log ( ^^^^) +2 log ( +loge/ - -8.81+2 log f ^ 



6GMe J \RE + hJ ■' \RE + h^ 

(53) 

As can be seen in Fig. (2), this leads to a new excluded region in the (a, A) parameter 
space. 

6.2.2. Periodic frequency shift 

As in the previous paragraph, one may obtain the periodic clock rate shift Eq. (49), 

.uGMe Re + h 2^2 /^.\ 
ep~-^^^2 ^^e^ = 2ece^ <^eo , (54) 

clearly below the constant contribution. 

6.2.3. Propagation time delay 

The additional propagation time delay reads 

GMEa h 

^ty- — • ^^^^ 

If the effect is undetected at a level At 10~^^ s, one obtains 

a<l^^^2.9xW-^(-^] . (56) 
GMe h V 1 my ^ ^ 

For a lower bound of A « 10* m (only one order of magnitude above RE,h), the 
result a < 0.1 is foimd, which does not advance the already excluded region of the 
parameter space (A, a) (see Fig. 2). 

As Fig. 2 shows, if clock rate shifts induced by a long-range Yukawa potential are 
undetected at the considered precision, then the excluded region already probed in 
the A = 10^ — 10* m domain becomes even deeper. This is a novel result stemming 
from this approach. 



7. Detection of a power-lav^f addition to the Newtonian potential 

One finally approaches the possibility of additions to the gravitational potential of 
the form 



.,.-^(fl . ,57) 
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where n 7^ —1 is a (possibly non- integer) exponent and i? is a characteristic length 
scale arising from the underlying physical theory. 

Phenomenologically, such a modification of the law of gravity is an interesting al- 
ternative to the more usual Yukawa parameterization, and allows one to investigate 
a wider range of extensions and modifications of GR. 

Such an addition can also be theoretically motivated: it arises from power-law 
induced effects at astrophysical scales due to the so-called Ungravity inspired sce- 
nario, which at short range involves the exchange of spin-2 unparticles of a putative 
scale invariant "hidden" sector within the Standard Model ^^'^o. Bounds for these 
"Ungravity" -inspired corrections can be obtained from stellar stability considera- 
tions cosmological nucleosynthesis and the gravitational quantum well . If a 
power-law addition is related to Ungravity, the exponent n follows from the scaling 
dimension of the unparticle operators du, through n = 2du — 2; the lengthscale R 
reflects the energy scale of the unparticle interactions, the mass of exchange particles 
and the type of propagator involved. 

Other possible power-law additions to the Newtonian potential may arise from 
f{R) theories of gravity 5^,55,56^ which generalize the Einstein-Hilbert action by 
considering a non-trivial scalar curvature term and/or a non-minimal coupling of 
geometry with matter ^^): the extra contribution arising from these models have 
been obtained in an astrophysical context, when addressing the puzzle of the galaxy 
rotation curves ^8,59^ .^gjj ^ -^^ ^ cosmological one, when applied to the issue of 
the accelerated expansion of the Universe 

From Eq. (57), one sees that the Newtonian potential <i?Ar is recovered by setting 
R = (for positive n) or i? — >■ 00 (for negative n). The limit n ^ is ill-defined, 
since the additional term Vp docs not vanish, but is instead equal to the Newtonian 
potential, Vp = ^n- hence, one should rewrite the gravitational constant in terms 
of an effective coupling, leading to the full potential 



where i?o signals the distance at which the full gravitational potential matches the 
Newtonian one, <l>(i?o) = ^n{Ro)- 

This additional length scale Rq should be an integration constant, obtained after 
solving the full field equations that lie behind the considered power-law correction. 
For simplicity, one assumes that (R/Rq)" <C 1, so that this term may be safely 
discarded — at the cost of neglecting the regime n ^ 0. With this in mind, the 
following paragraphs use Gp = G freely; one remarks that this approach is comple- 
mentary to that considered in Ref. 51. 




(58) 



with 




G 



(59) 
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7.1. Constant frequency shift 



As before, the constant relative frequency shift Cc of an emitted signal is given by 
Eq. (21), which is dominated by the linear Doppler shift, 

where the dimensionless ratio ^ = R/Re is defined. Comparison with a precision 
e/ 10~^^ yields the constraint 

which gives an upper or lower bound on ^, depending on the sign of n. Since the 
l.h.s. grows with n, one concludes that ^ ^ 1 + H/Re ~ 4.64 for large n. The 
allowed region in the parameter space (n, R) is depicted in Fig. 3. 

Analogously to the discussion of the previous section, the range e < 1 ^ R < Re 
would imply the inclusion of a suitable form factor to model the exchange of spin-2 
unparticles within the extended mass distribution of the Earth. The results of Fig. 
3 are obtained with no form factor, thus presenting an unrealistic upper bound for 
the (n, ^) parameter space: the inclusion of a suppressing form factor would only 
increase the allowed region for positive n and ^ < 1, which is mostly unrestricted. 

7.2. Propagation time delay 

As before, one also computes the additional propagational time delay, 



At, = 



( Re 



1.48 X 10""— (0.22" - 1) s . (62) 

n 



\RE + h^ 

This expression may be simplified by considering the large positive or negative n 
regimes. 

7.2.1. Large positive n 

For a sufficiently large positive exponent n, Eq. (62) reads 

Atp ~ ~ -1-48 X 10""— s , (63) 

and comparison with the time resolution At ~ 10~^^ s yields the upper bound 

— < 6.76 X 10"^ . (64) 
n 

If ^ > 1, the r.h.s. grows and this inequality is quickly violated, similarly to the 
previous paragraph. Hence, one obtains the upper bound ^ < 1 for n 0, as 
depicted in Fig. (3). 



July 19, 2011 0:17 WSPC/INSTRUCTION FILE GalileoIJMPD4 



Using global positioning systems to test extensions of General Relativity 23 

7.2.2. Large negative n 

If n <C 0, Eq. (62) simplifies to 

A^P - ^ 1.48 X io-"i£:^ . . (65) 

nc^ \RE + hJ n ^ ' 

Equating this to the considered time accuracy 10~^^ s of Galileo, one obtains the 
inequality 

i^:^ « 6.76 X 10-2 , (66) 
\n\ 

As before, the l.h.s. grows with n if 0.22^ > 1; thus, one recovers the previously 
obtained bound f > (0.22)"^ ~ 4.64 for large n. Notice that this result is indepen- 
dent of the time resolution, since a large negative n regime eventually fulfills this 
condition, for any given Af. The allowed region for the ^, n parameters is depicted 
in Fig. 3, and this asymptotic bound is visible in the negative n axis. 

8. Conclusions 

In this work, one has assessed the possibility of detecting signals of new physics 
through the use of the Galileo navigation system. This application could be valu- 
able, as any unexpected new phenomenology could provide further insight into what 
lies beyond the Standard Model of particle interactions and GR. One has specifi- 
cally looked at the propagation time delay and frequency shift induced by several 
corrections to the Newtonian potential: PPN second-order terms, a cosmological 
constant, a constant anomalous acceleration, a Yukawa addition and a power-law 
addition. 

As it turns out, there is no possibility of using the Galileo satellite naviga- 
tion system to currently detect PPN second-order additions nor a cosmological 
constant with the reported value; likewise, a constant acceleration of the order 
10"^ — lO"^*^ m/s^, characteristic of the Pioneer anomaly or MOND, also falls well 
below dctectability. 

One finds that the excluded region for the parameter space (A, a) of a Yukawa 
interaction can be extended with the available precision: the A ~ 10^ — 10^ region 
is deepened so the allowed range for a is lowered by about one order of magnitude, 
to a < 10-8. 

Moreover, Fig. 2 shows that a possible future generation of satellites, possessing 
onboard clocks with a relative frequency accuracy 10"^^, might probe a relevant, 
yet unassessed zone of this parameter space; this kind of clock stability will be 
eventually available with the onset of quantum time synchronization and subsequent 
space qualification. 

Another interesting result of this study lies in the exclusion of a definite range 
of values for the characteristic lengthscale of a power-law addition, R, namely those 
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Fig. 3. Allowed regions (gray) for a Ungravity-type additional force with lengtliscalc R and 
exponent n, and superimposed limit obtained for a relative frequency precision e/ = 10^ (full) 
or e = 10~^® (short dash), and time resolution At = 10~^^ s (long dash) or At = 10~^^ s (dotted). 



lying between < i? < 4.64Re- For positive exponents n, the region R < Re is 
essentially unconstrained (with the allowed region further enlarged by the unac- 
counted suppression due to the inclusion of a form factor). 

Despite of being somewhat incomplete and exploratory, this study serves to show 
that the Galileo navigation satellite system could prove to be a valuable instrument 
for improving our understanding of fundamental physics. This result is encouraging, 
and can be considerably improved in a future generation of Galileo — particularly 
if well-thought science objectives are incorporated into the design of the ground 
segment and satellite constellation. 

Furthermore, results concerning the periodic clock rate shift arising due to de- 
viations from a perfectly circular orbit show that a dedicated satellite with a state 
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of the art onboard clock and a highly eccentric orbit could serve to probe yet un- 
constrained domains of the parameter space of the discussed extensions of GR. 

Clearly, this goal is not within the framework set by any global navigation satel- 
lite system, which by design aim at a negligible eccentricy orbit. However, one can 
envisage that a small subset of the constellations could be launched with a suffi- 
ciently large amount of fuel (adding to the already provided station-keeping and 
controlled orbital decay supply), so that its eventual deorbiting could first involve 
some sort of elliptical path: such dual-use could provide interesting scientific oppor- 
tunities in a cost-effective way (when compared with a dedicated probe), and serve 
to highlight the global value of a particular system, beyond its basic navigation 
purpose. 

Although not within the scope of this study, the Galileo satellite navigation 
system could also prove to be highly valuable in probing the validity of the Local 
Positioning Invariance (LPI) principle This postulate, one; of the fundamental 
pillars of General Relativity, asserts that clock rates are independent of their space- 
time positions; the available experimental constraints tell us that this invariance 
holds down to a relative level of 1.4 x 10~^ ^^'^^ (see Ref. 63 for the most recent 
bounds, and Ref. 64 for a discussion). 

Endowing one or more elements of the Galileo constellation with more precise 
clocks and providing sufficiently stable communications with ground stations (pos- 
sibly through a microwave link) , an improvement of up to two orders of magnitude 
on the LPI could be achieved. However, this could be somewhat difhcult, given that 
one does not have direct knowledge of the emitted frequencies, only those of the 
received signals and the master clocks on the ground. 

Indeed, the frequency shift inferred from the comparison between these and the 
broadcasted timing information could suffice (with increased precision) to ascertain 
the modifications of gravity here addressed, because one can use the multiple satel- 
lites available to correct for some unknown instabilities in the emitters, in a sense 
obtaining an "averaged" signal. However, measuring the LPI would involve compar- 
ing the individual space-borne clocks and looking for deviations inferred from the 
received signals, and it could prove difficult to distinguish actual clock-rate differ- 
ences due to a violation of the LPI from simple drifts from the original base rate 
(which should nevertheless be small, for the 10~^^ stability). 

To counteract this issue, an alternative could involve the use of cornercubes on 
the surface of one or more elements of the Galileo system, thus enabling the use of 
accurate laser ranging: the obtained redundancy would help to lift this degeneracy 
and signal "true" violations of the LPI. Such an extension of the currently envi- 
sioned system could be dubbed Siderius Nuncius, the Celestial Messenger, after the 
1610's historic account of Galileo of the first use of the telescope for astronomical 
observations. Siderius Nuncius might constitute a valuable instrument to advance 
our understanding of the mysteries of the Cosmos. 
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